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OBJECTIVE 


The objective of this discussion is to examine the application of 
eigenspace design techniques to an active flutter suppression system for 
the DAST ARW-2 research drone. Eigenspace design techniques allow the control 
system designer to determine feedback gains which place controllable eigenvalues 
in specified configurations and which shape eigenvectors to achieve desired 
dynamic response. Eigenspace techniques have been applied to the control of lateral 
and longitudinal dynamic response of aircraft [1,2]. However, little has been 
published on the application of eigenspace techniques to aeroelastic control 
problems . 

This discussion will focus primarily on methodology for design of full- 
state and limited-state (output) feedback controllers. We do not intend to address 
the significant difficulties associated with the realization of full- 
and limited-state controllers. Most of the states in aeroelastic control 
problems are not directly measurable, and some type of dynamic compensator is 
necessary to convert sensor outputs to control inputs. Compensator design can 
be accomplished by use of a Kalman filter modified if necessary by the Doyle- 
Stein procedure for full-state loop transfer function recovery [3], by some 
other type of observer, or by transfer function matching. 
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EIGENSPACE DESIGN TECHNIQUES 


Eigenspace techniques allow the designer to place closed-loop eigenvalues 
(A^) and shape closed-loop eigenvectors (v^) . We will briefly review the theory. 
For a more detailed discussion see Refs. 1, 2, and 4. First we assume the system 
is controllable and observable and the matrices B and C are of full rank. (In 
the case of full-state feedback, C = I.) Later the controllability assumption 
will be relaxed. The above assumptions yield the results shown below. If we 
have full-state feedback and n controls, we can arbitrarily place all eigen- 
values and shape all eigenvectors to any desired form. If we have full-state 
feedback and a single control, only pole placement is possible. Since any 
attainable eigenvector is in the subspace spanned by (A I-A)“1 b, it is 
impossible to exactly achieve a desired eigenvector in most aircraft control 
problems. In practice this does not appear to be a serious problem. 

Uncontrollable eigenvalues cannot be moved but an additional element in each 
eigenvector associated with these eigenvalues can be shaped. 


x = Ax + Bu 
y = Cx 
Dim (x) = n 
Dim (u) = m 
Dim (y) = r 

(1) max (m,r) closed-loop eigenvalues can be assigned 

(2) max (m,r) closed-loop eigenvectors can be shaped 

(3) min (m,r) elements of each eigenvector can be arbitrarily chosen 
Attainable eigenvector in space spanned by 

(IA.-A) _1 B 
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CALCULATION OF GAIN MATRIX-I 


We will first describe the eigenspace design technique for full-state feed- 
back for a system described in standard state-space form. The design procedure 
consists of determining a gain matrix K such that for all closed-loop eigen- 
value and eigenvector pairs 

(A+BK)v^ = A^v^ 

where A^ is the desired closed-loop eigenvector and v^ is the associated closed- 
loop eigenvector. This is equivalent to finding w^ such that 

(IA^-A) v ± = Bw i 

Once all the w^'s have been found, the gain matrix can be calculated. 

In order to arbitrarily place all the A.'s and v^'s, the control vector will 
have to be of the same order as the state vector and B would have to be inverti- 
ble. In general this is not the case, and the achievable eigenvalues must lie 
in the subspace spanned by 

(A i I-A)~ 1 B 

In general, the desired eigenvector v^ will not reside in this subspace. 


x = Ax + Bu 


y = Cx 


u = Ky 

x = S a v.e'V 
i-1 1 1 


(A + BK)v.=A v 
1 i l 


(IA . -A)v .=Bw . 
x xx 


v 

i 


L .w- 
x x 


L. = (A.I-A) 1 B 
x x 

K = WCCV] -1 
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CALCULATION OF GAIN MATRIX-II 


Since the desired eigenvalues are in general not achievable, the w.’s are 
selected to minimize the weighted sum of the squares of the difference between 
the elements of the desired and attainable eigenvectors given by the perform- 
ance index J-^. The term P^ is a positive definite symmetric matrix whose 
elements can be chosen to weight the difference between certain elements of 
the desired and attainable eigenvalues more heavily than others. Setting the 
derivative of with respect to w^ equal to zero gives w-^. The notation * 
denotes complex transpose. Once w^ is calculated, the achievable eigenvector 
v^ is obtained. If an eigensolution is not to be altered, setting w-j^ = 0 
assures that the associated v^ and remain in their open-loop configurations. 
If we desire output feedback, the procedure is easily modified as shown. 

In case of a complex eigenvalue, a real-gain matrix results from a simple 
transformation [1, 2, 4]. 


J. 

i 



* d 

- v.) P.(v. 

l l i 


V 


3J. 
i 

3w. 

i 


0 


* —1 * H 

w. = (L. P.L.) L. P.v. 

l ill ill 


If an eigensolution is not to be altered, 
w. = 0 
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UNCONTROLLABLE EIGENVALUES-I 

In aeroelastic control problems, the states associated with the gust model 
are uncontrollable. Moore [4] showed that it is possible to use feedback to 
assign some components of eigenvectors associated with uncontrollable eigen- 
values. An algorithm for performing such an assignment is given. For 
an uncontrollable eigenvalue A , 

[I A -A] 

8 

is singular. We can partition the eigenvector Vg associated with this 
eigenvalue as shown below where v E ^ contains only the uncontrollable states. 


A uncontrollable 
g 


Partition v such that 
g 


[A 

g 


I-A]v =Bw 
g g 


Becomes 



v contains only uncontrollable states 
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UNCONTROLLABLE EIGENVALUES -II 


The equation 


Rv 


II 


0 


o IX 

is automatically satisfied if we select v g to be equal to the open-loop 
portion of v g which contains the uncontrollable states. Since A g is not an 
eigenvalue of A 1 , [A I -A 1 ] is nonsingular, and by performing the indicated 
calculations, w g can g be determined in much the same way as for the controllable 
eigenvalues. 


Rv =0 

g 

T T -1 T I -1 II 

v = [A I-A ] B w - [A I-A ] Pv 
g g g g g 

v ^ = L w + Av 
g g g g 


minimizing 


. Id I-* , Id In 

J = (v - v ) P (v - V ) 
g g g g g g 


yields 


w = (L *P L ) 1 L *P (v * d - Av ) 
g ggg ggg g 
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DAST ARW-2 SENSOR AND C 


SENSORS 


/ERTICAL ACCELEROMETERS 
WBL 84 REAR SPAR 
WBL 82 FRONT SPAR 
: WBL 92 REAR SPAR :i 
WBL 92 FRONT SPAR. 

• FSS 


VERTICAL ACCELEROMETER - 
BS 250 



• GLA 
t AFCS 


OARD AILERON 


ILIZER 


ROLL, PITCH, YAW RATE - 
ANGLE OF ATTACK (NASA) 

• RSS 

• AFCS 


D AILERON 








PERFORMANCE SPECIFICATIONS 


The design flight condition for the flutter control system is M=0.86 
(275 m/s or 908 ft/s) and an altitude of 15000 ft (475 m) . At this flight 
condition, the uncontrolled wing flutters, and the flutter control system is 
required to stabilize the wing without exceeding specified limits on rms control 
surface activity. The control surfaces saturate if these limits are exceeded. 

Gain and phase margins must be adequate. The wing flutters at M=0.75 at this 
altitude, and the flutter control system must be activated at M=0.7. It must 
be verified that activation of the control system does not destabilize the wing 
at this flight condition. Also the flutter controller must not result in excessive 
increases in bending, shear, or torsional loads compared with the uncontrolled wing. 


Design Condition 

M = 0.86 h = 15000 ft. 


Maximum RMS Control Surface Activity for 12 ft/s Gust 



Deflection 

Deflection Rate 

Inboard 

10° 

130°/s 

Outboard 

15° 

740°/s 


Mininum Stability Margins 
Gain - 6 dB 
Phase - 45° 

No large increases in bending, torsion, and shear at M = 0.7 


249 



AEROELASTIC MODEL 


In the aeroelastic model of the wing given below, y^ is the vector of 

displacements of the various flexural modes, y is the vector of control surface 

deflections, y is the gust velocity, M is the structural mass matrix, C is the 
g s s 

structural damping matrix, is the structural stiffness matrix, q is dynamic 

pressure, and Q c is the matrix of aerodynamic influence coefficients. Q c is 

calculated as a function of reduced frequency by a doublet-lattice procedure 

and is approximated by Q^, a matrix of rational polynomials in the Laplace operator 

s. The matrices A^ are selected to give the best least-squares fit to Q c over a 

range of reduced frequencies. 


( CM ]s 2 + [C ]s + [K ] ) [y ] + q [Q (s) ] 
s s S t c 


[Q c Cs)]-tA 0 ] + [A-^] 2v + 

c m 


cs 

2v 


J_i 

+ £ 

7T1=1 


rSp+2] s 
s+ 2v v 
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FLEXURAL MODES 


Initially, seven structural modes were used in the math model of the wing; 
however, by comparing eigenvalues calculated using lower order models with 
eigenvalues resulting from a model which included seven structural modes, it was 
found that flutter could accurately be modeled by including only three modes. 
These modes are. shown below. The first mode will be labeled first mode bending, 
the second mode will be labeled second mode bending although it contains some 
torsion, and the third mode will be labeled first mode torsion. Rigid-body modes 
were not included . 




Second Bending 



251 



I 


OPEN-LOOP ROOT LOCUS 


The locus of the roots associated with the flexure modes is shown below. 

The lowest frequency mode is associated with first mode bending, the next highest 
with second mode bending and the highest with first mode torsion. It can be seen 
that first mode bending is the unstable mode while second mode damping increases 
with velocity. The frequencies of these modes approach one another with increas- 
ing velocity. The first torsion mode is not affected by velocity as much as the 
other two modes, but it is necessary to include this mode in order for the system 
to flutter. The wing flutters at a velocity of 787 ft/s (240 m/s). This is about 
Mach 0.75. 
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WIND GUST AND ACTUATOR MODELS 
AND STATE SPACE FORMULATION 


A second-order model forced by white noise was used to simulate the vertical 
gust. Both inboard and outboard ailerons are driven by high bandwidth actuators. 
In the range of frequencies covered by the three-mode structural model, a fourth- 
order transfer function was shown to give a very close approximation of the actual 
inboard actuator/aileron transfer function. A third-order transfer function was 
used for the outboard aileron. The details of these models are given in Ref. 5. 
The state space model of the combined system is given below. The vector Xp 

includes the displacements and velocities associated with flexure modes and the 
aerodynamic lag states. The vector X c includes the states associated with the in- 
board and outboard actuator models. The vector Xg includes the states associated 
with the gust model. The vector U is the control input to the actuators, and w is 
the scalar white noise input to the gust model. The total system model is 18th 
order. Note that the open-loop responses of the actuators are decoupled from one 
another and are not influenced by the motion of the wing (small inertial cross- 
coupling terms have been neglected) and that the gust states are uncontrollable. 


V" ' G g (s) 


y c /'-Ui = C (s) 
i 


y c /u o ~ G o (s) 




A 11 

A 12 

A 13 
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0 
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X 
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+ 

B 

U + 
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A 33 
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EIGENSPACE DESIGN APPROACHES 


The initial eigenspace controller was designed by rotating the unstable eigen- 
values about the imaginary axis and leaving all other eigenvalues and eigenvectors in 
their open-loop positions. The results are shown in the table on the following page. 
Although this initial design stabilized the wing at both the flutter test condition 
(M=0.86) and at the condition at which the flutter controller would be initially 
activated (M=0.7), the rms inboard deflection rate is near its maximum (saturation) 
value at the flutter condition. It was felt that the performance might be enhanced 
by redesign of the control system to reduce the inboard deflection rate. Also the 
initial design approach did not use the capability of eigenspace techniques to shape 
eigenvectors as well as assign eigenvalues, and it was desired to exercise this 
capability. Since the aircraft exhibits satisfactory response at velocities less 
than the flutter speed, it was decided to force the closed-loop response of the 
wing at the design condition to approach the open-loop response of the wing at a 
velocity 656 ft/s (20% less than the flutter speed). The closed-loop eigenvalues 
associated with the flexure modes and aerodynamic lag states were moved to their 
open-loop positions at 656 ft/s. The open-loop actuator eigenvalues were the same 
for both flight conditions and were not moved, and the gust eigenvalues were uncon- 
trollable and could not be moved. The desired eigenvectors were selected to be the 
open-loop eigenvectors of the wing at 656 ft/s. The weighting matrices in the per- 
formance index were initially set at one. This procedure reduced all surf ace activ- 
ity only slightly. It was decided to use eigenvector shaping to shift control 
surface activity from the inboard to the outboard actuator. The components of the 
open-loop aeroelastic eigenvectors in the actuator directions are zero for all 
flight conditions. Since the desired eigenvectors were chosen to be the open eigen- 
vectors, penalizing the difference between the achievable and desired aeroelastic 
eigenvectors in the direction of the inboard actuator would reduce inboard activity. 
This was accomplished by increasing the weights on these components to 2.5 x 10^ 
while all other weights remained at one. The results are shown in the table on the 
following page. 
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COMPARISON OF RMS CONTROL SURFACE ACTIVITY FOR 
THE EIGENSPACE CONTROLLERS 


The performance of the initial and final eigenspace controllers is summarized 
in the table below. The inboard actuator rate for the initial eigenspace design 
is almost saturated while, with eigenvector shaping, the control effort is shifted 
from the inboard to the outboard control surface. Failure of both the inboard 
and outboard aileron was simulated. When the Inboard actuator failed, performance 
was not affected very much; however, failure of the outboard actuator resulted in 
an unstable response. This indicates that the inboard actuator is not a good 
choice for use in flutter suppression. 


Controller Design 

* 

Inbd Defl 

Inbd Rate 

Outbd Defl 

Outbd Rate 






Unstable 
Roots Rotated 
About Imag Axis 
(Initial Design) 

0.5 

108.0 

3.3 

509.0 

Eigenvector Shaping 
to Reduce Inbd Rate 
(Final Design) 

0.9 

86.0 

4.7 

612.0 

Initial Design 
Inbd Failed 


- 

3.5 

519.0 

Final Design 
Inbd Failed 

- 

- 


572.0 

Max Allowable 

10 

130 

15 

740 
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CLOSED-LOOP ROOT LOCUS 


The root locus of the aeroelastic modes for the final eigenspace design is 
shown below. The wing goes unstable at about 1017 ft/s (M=0. 96) . The open-loop 
flutter speed is 787 ft/s (M=0.75); therefore, the control system results in an 
increase in flutter speed of about 29%. As in the open-loop case, the first 
bending mode goes unstable, but in the closed-loop case, the eigenvalues associated 
with this mode move to the real axis where one real root goes unstable. The 
roots associated with the second bending mode are almost unstable at this velocity. 
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FLUTTER BOUNDARY 


The results of varying altitude while maintaining Mach number constant are 
used to define the flutter boundary for the open-loop wing and the wing controlled 
with the final eigenspace design. At M=0.86, the uncontrolled wing is unstable 
until an altitude of 6700 m is reached. At the same Mach number, the controlled 
wing is stable for altitudes above 2900 m. At M=0.7, the uncontrolled wing is 
stable for altitudes above 1800 m, whereas the controlled wing is stable for 
altitudes above 2100 m. 



257 



STABILITY MARGINS 


Since the inboard aileron was ineffective in stabilizing the wing, and the 
outboard aileron is critical, stability margins with only the outboard loop closed 
are shown below. The initial eigenspace design results from rotating the unstable 
eigenvalues about the imaginary axis and, with the inboard aileron inactive, this is 
identical to a design resulting from linear quadratic regulator theory. This is 
guaranteed to have excellent stability margins (Ref. 6). The gain margin is 6 dB, 
and phase margins are greater than 60° . The frequency response characteristics of 
the final eigenspace design are considerably different from those of the initial 
design. Gain margins are 6 dB or better, but phase margins are less than 20°. The 
solid lines represent the initial eigenspace design, and the dotted lines represent 
the final eigenspace design. 
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OUTPUT FEEDBACK 


Since the inboard actuator was ineffective for flutter control, it was elimi- 
nated from the design. This resulted in a system with a single control input and 
allowed only eigenvalue assignment. An accelerometer was used to measure the 
motion of the wing, and a feedback compensator that approximated the frequency 
response characteristics of the full-state loop transfer function was designed. 
Also the effects of restricting the number of states in the feedback controller 
were investigated. The design approach was to treat this as a problem in output 
feedback as described earlier. The C matrix was selected to eliminate various 
states from the output, and a feedback controller was designed using eigenvalue 
placement techniques. Since the gust states are uncontrollable if P states are 
fed back, only P-2 eigenvalues could be placed. The following sets of states were 
eliminated from being fed back: (a) gust (2 states), (b) aerodynamic lags 

(3 states), (c) actuator (3 states), (d) first bending (2 states), (e) second 
bending (2 states), (f) first torsion (2 states), and (g) first torsion and 
aerodynamic lags (5 states) . Except in the case where the first bending mode 
was eliminated, unstable eigenvalues were rotated about the imaginary axis, 
and eigenvalues associated with the retained states were maintained in their 
open-loop positions. The positions of the other eigenvalues were not assigned. 
Since the eigenvalues of the first bending mode were unstable, these eigenvalues 
were rotated, and the eigenvalues of the first torsional mode were not assigned. 
This mode then went unstable. The results are summarized below. The rms response 
was very sensitive to the gust states, and attempts to remove these states or 
alter their eigenvectors resulted in large rms responses. Loads at the wing root 
were reduced using the full-state design. Note a slightly different gust model 
was used so the rms results below are not directly comparable with those 
presented earlier. 


States 

Eliminated 

Otbd 
Def 1 

Otbd 

Rate 

Gain 

i 

Margin Db 

Phase Margin 
(Degrees) 

Gust 

32 

456 

6.0 

±60 

Actuators 

2.0 

265 

4.1 

±36 

Aero. Lags 

6.6 

394 

5.7 

-50,60 

Mode 1 

Unstable 

Unstable 

Unstable 

Unstable 

Mode 2 

2.0 

260 

2.8 

-27,18 

Mode 3 

2.0 

271 

5.8 

-63, +54 

Mode 3 + 
Aero. Lags 

6.6 

430.0 

6.1 

-54, +72 
+81, +45 

Full State 

2.0 

259 

6 

±60 

Compensator 

2.72 

250.8 

6.0 

±60 


RMS Responses M = .86 
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CONCLUSIONS 


Eigenspace techniques can provide a powerful tool for the design of feed- 
back control systems for aircraft. The design techniques described in this 
paper are easily implemented and are computationally inexpensive. The basic 
problems facing the designer are those of determining where to place eigenvalues 
and of selecting appropriate eigenvectors. This usually requires some insight 
into the system to be controlled; however, in this respect all control system 
design techniques are the same. Stability margins must be carefully examined 
since there are no guaranteed margins as is the case for linear quadratic 
regulators. On the other hand, modal decoupling is easily achieved, and certain 
roots can be maintained in their open- loop configurations if desired. Output 
or limited-state controllers can easily be designed. 


! 

Shear 

Torsion 

Bending 


(lbs x 10 2 ) 

(in - lbs x 10 2 ) 

(in - lbs x 10**) 

Open Loop 

4.131 

6.928 

2.365 

Closed Loop 

3.775 

1.029 

2.321 


Wing Root Loads M = 0.7 
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